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We present a modification of tfie generalized Flory dimer theory to investigate the nematic 
(N) to isotropic (I) phase transition in chain fluids. We focus on rigid linear fused hard-sphere 
(LFHS) chain molecules in this study. A generalized density functional theory is developed, which 
involves an angular weighting of the dimer reference fluid as suggested by decoupling theory, to 
accommodate nematic ordering in the system. A key ingredient of this theory is the calculation 
of the exact excluded volume for a pair of molecules in an arbitrary relative orientation, which 
extends the recent work by Williamson and Jackson for linear tangent hard-sphere chain molecules 
to the case of linear fused hard-sphere chains with arbitrary intramolecular bondlength. The present 
results for the N-I transition are compared with previous theories and with computer simulations. 
In comparison with previous studies, the results show much better agreement with simulations for 
both the coexistence densities and the nematic order parameter at the transition. 



I. INTRODUCTION 



The hard-sphere chain fluid has been widely studied in recent years as a model for polymeric and related liquids. 
Various theories have been applied to determine the thermodynamic properties of the model fluid, such as thermo- 
dynamic perturbation theory (TPT) generalized Flory dimer theory (GFD) |^-||, scaled particle theory (SPT) 
P-pl| and various integral-equation theories |12 1^]. Despite these studies, several limitations still characterize most 



theories for the hard-sphere chain fluid, including the independence of intramolecular structure on the surrounding 
medium in homogeneous fluids [§|j3j, the inability to account for liquid-crystalline ordering |||,^,^,^ and a general 
difficulty in formulating practical extensions of these theories to non-uniform systems. Several density-functional 
approaches |16|-p^ to the latter problem have been examined, but these have not overcome the other limitations 
mentioned above. 

In this article, we investigate a systematic and straightforward density-functional method to begin resolving these 
limitations. We focus on a rather specific problem, namely to characterize the nematic-isotropic (N-I) phase transition 
in a fluid composed of rigid linear fused hard-sphere (LFHS) chain molecules. Several related works have examined 
this topic recently ||,|ll|,^,^-|2l[l . Mehta and Honnell (|] applied GFD theory to the LFHS fluid, but were restricted 
to examining the isotropic phase. This is also true of the integral-equation and SPT approaches studied in ref. [^ , 
although some information about the N-I transition could be deduced from stability analysis of the isotropic phase. 
Refs. [0,|9|-|2^ studied the N-I transition in linear tangent hard-sphere (LTHS) chain fluids, using a density-functional 



approach based on the decoupling approximation introduced by Parsons and Lee 1 23 1 for hard ellipsoids and hard 
spherocylinders. This approximation consists of a rescaling of the virial series originally used by Onsager to 
account for the thermodynamic properties of a rigid-rod fluid. (Ref. ||2l| also describes an extension of the method to 
a fluid of semi-flexible chains.) 

The formulation discussed here is very similar to that of refs. but utilizes a concept closely related to that 

of the TPT and GFD theories, namely that the thermodynamic properties of the full chain fluid can be described 
by a judicious combination of those characterizing "reference" fluids composed of smaller sub-units, in particular 
monomers and dimers. Both formally and intuitively, this picture is well suited to describing fluids of flexible chain 
molecules [|]-|7|. A key ingredient of the present work, focusing on rigid rods, is that the excluded volume (and hence 
the second virial coefficient) for a pair of LFHS molecules in an arbitrary orientation can be written exactly as a linear 
combination of that for a pair of diatomic hard-sphere molecules and a pair of hard-sphere monomers. This property is 
based on a recent analysis of the orientation-dependent excluded volume for LTHS chains by Williamson and Jackson 
p5| . Here it is shown that the results of ref. |^ can be extended straightforwardly to linear fused hard-sphere chains 
of arbitrary intramolecular bondlength. 

In the present article, the decoupling approximation is utilized in two stages. In the first stage, it is assumed that all 
higher-order virial coefficients scale (in terms of monomer and diatomic contributions) in the same way as the second 
virial coefficient. In the case of an isotropic fluid phase, the resulting approximation for thermodynamic properties 
is similar to that of GFD theory, but differs from the latter in reproducing the exact second virial coefficient. In the 
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second stage, the excess free energy of the reference dimer fluid is rescaled in terms of the free energy of the isotropic 
phase of this fluid. This rescahng incorporates the exact weighted angic-avcrage of the pair excluded volume, as in 
the theories of refs. ||l9|~[2l||, and thus is able to account for nematic ordering, in contrast to the original GFD theory. 
The theory is p resented in Section ||. Although the focus of this paper is on uniform phases of a rigid molecular 



fluid. Section [I A begins with the more general case of a non-uniform fluid and extensions of the theory to semi-flexible 



molecules are indicated briefly. In Sections [IB and [IC, the theory is specialized to uniform phases of LFHS chains. 
Calculations and results are presented in Section III, including comparisons with available Monte Carlo data for the 
N-I transition of both LTHS and LFHS chains. It is found that the present theory overestimates the reduced pressure 
when compared with both simulations and previous theories. However, the present results for coexistence densities 
and the nematic order parameter at the transition are in excellen t ag reement with simulations and arc more accurate 



than those predicted by other recent theories [|liyiSH2l[|. Section IV contains a summary and some conclusions 



II. THEORY 
A. General 

Here we consider a one-component fluid containing on average N molecules, each molecule (or "n-mer") consisting 
of a rigid linear chain of n atomic sites. The configuration of any such molecule labelled i is specified by the position 
Ti of some arbitrarily chosen "center" of the molecule and by the orientational Euler angles {di,4>i) = uji of the 
molecular chain axis relative to some space-fixed frame. We assume that the total intermolecular potential energy for 
N molecules is pairwise additive, and denote the pair potential between molecules i and j as C/(2) (r^Wi, rjWj). 

The probability density for finding any molecule in a configuration ruo is denoted p(ra;), which is normalized such 
that 

j drdw p{yuj) = N. (1) 

The integration over molecular position is constrained to occur within the system volume V . In density- functional 
theory, the Helmholtz free energy F of the fluid is expressed as a functional of p{vLd). A formally exact, generalized 
virial expansion of F in powers of p{rLo) can be obtained by straightforward adaptation of that for monatomic fluids 
p6[ , and is given by 

F ^kT j dYdujp{YLo)[\n{ATTp{ruj)K) - 1 + pUexivLo)] -f AF, (2) 

where T is temperature, k is Boltzmann's constant, (3 ~ (kT)^^, A is the thermal de Broglie "volume" of a molecule, 
and Uex is any external potential fleld acting on the fluid. AF is the excess (over an ideal gas) Helmholtz free energy 
of the system, and is expressed by the virial series 

oo 

AF = kTy (3) 

m—2 

The generalized virial coefficients _B,„ have a standard diagrammatic representation in terms of irreducible cluster 
integrals, where the vertices of the diagrams represent the product of functions /^(riwi), p(r2ti^2), ■ • • p(rmt^m) and the 
bonds denote the Mayer function 

/(i,j) = e-''^(^)('J)-l. (4) 
For simplicity, the label i is used to denote ViLJi. Thus, 

B2 = ~\l dld2p{l)p{2)f{l,2), (5) 
^3^~IJ dld2d3pil)p{2)pi3)fil, 2)/(2, 3)/(3, 1), (6) 

etc. 

The preceding equations can be formally generalized to fluids of semi-flexible molecules on interpreting the symbol 
i to include all internal coordinates of molecule i and replacing Uex{i) by the full one-body potential incorporating 
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internal bonding and flexibility constraints ^7|j28| . In the usual approach of density-functional theory, the distribution 
function p(i) in the preceding expressions for F is considered to be an arbitrary function of the molecular configuration 
i, while the equilibrium form of that distribution function is obtained by functional minimization of the corresponding 
grand canonical potential = F — Nfi, where fi is the chemical potential. 

In the standard application of the decoupling approximation for the case of uniform fluids ^2 23 2^ , the 2nd virial 
coefficient is treated exactly while the mth-order virial coefficient is approximated as 

= , (7) 

^2 

where is the virial coefficient of an appropriate reference fluid. In previous works, the reference fluid has been 
chosen to be either a hard-sphere fluid [^,^ or the isotropic phase of the actual system being considered p9| , p0| , p9[ . 
In either case, resummation of the virial series d) then gives 

AF = ^^AF'^^'f. (8) 

(A generalization of this procedure for a non-uniform fluid of hard spherocylinders is described in ref. |30[| .) The 
GFD and TPT theories for fluids of chain molecules are based on a different choice of reference system, involving the 
decomposition of each molecule into its constituent monomers and constituent dimers. Although originally derived 
using quite different arguments, we can obtain the GFD theory using the present framework. First, at the level of the 
2nd virial coefficient, we postulate that B2 can be written as 

S. = ar)sr+a(^)sf), (9) 

where al*^'* and a^''^ are appropriate mixing parameters. Here, b'',^^^ is the generalized 2nd virial coefhcient charac- 
terizing the fluid which results on decomposing each chain molecule into its constituent monomers, while -62^'' is the 
analogous quantity for the reference fluid of constituent dimers. Formally, the latter decomposition assumes that the 
number of monomers per chain n is even. However, we anticipate (as suggested by previous work on GFD and re- 
lated theories P,[7|,^ ) that both the monomer and dimer reference fluids should be characterized by suitable effective 
values of the number of monomers um and number of dimers no , respectively. In accordance with the decoupling 
approximation, we assume that the mixing parameters al*^'' and can be chosen so that (p|) generates the exact 



B2 ■ This condition will be examined in subsection [IB. If it is assumed that all higher-order virial coefficients B^ 
scale in a similar manner, i.e. 

B„ = a(*^)s(?'^)-fa(^)s(f), (10) 

then resummation of the series in (^) gives 

AF = ai*^)AF(*^)[pM(r)] + a^^) AF^^) [pz5(rc.)]. (11) 

In this equation, pM(r) and pd{vuj) are the probability densities for monomers and dimers, respectively, while AF*^*^-* 
and AF(^) are the corresponding excess free energies, which in general are functionals of the densities. For simplicity, 
we assume that the n-mer contains only one type of monomer and dimer sub-unit; the generalization to heteronuclear 
molecules is straightforward. 

The standard GFD theory for isotropic, homonuclear hard-sphere chain fluids, as investigated previously in refs. 
|^-||, can be formulated from the theory presented here. In particular, GFD theory is obtained from eqs.(^), (pl|), 

and a modified form of the condition (|^). In order to determine the parameters al*^^ and a^^^ in this case, ^ is 

applied with B2 and B^^ approximated by the cross second virial coefficient between a monomer and the full 71-mer 

and between a monomer and a dimer, respectively. An additional condition is then required to fix ai^^'' and al^'' 
uniquely, which turns out to be 

nMaW+nz,al^) = l. (12) 

This is consistent with the "initial conditions" that ai*^'* — 1, a^ri'^ = when n = um = 1 and al*^'' = 0, ai^-* = 1 
when n — um = 2n]j = 2. With these approximations, one obtains 



riMa. 



L*''=-r„ , nDai^) = l + F„, (13) 



3 



where Yn is the quantity defined in eq.(14) of Honncll and Hall H]. 

In the next subsection, we shall specialize and ( |ll| ) to the case of a uniform fluid of rigid homonuclear hard- 
sphere chains. It will be shown that (0) can be satisfied exactly when using the exact values of B2 and i?2^'' and with 
a unique choice for the parameters a„ •* and ai^^ depending only on the geometric properties of the molecules. 



B. Application to LFHS Chains 



Here we consider that the full chain molecule is a rigid rod composed of n identical hard-sphere atoms, each of 
diameter d, with adjacent atoms separated by bondlength I. In a uniform but possibly orientationally ordered fluid, 
the n-mer probability density takes the form 

(14) 

where p = N/V is the molecular number density and /(w) is a normalized angular distribution function. The 
generalized second virial coefficient B2 defined by ^ becomes 

p^V f 

B2 = — j dri2dLJidLJ2f{uJi)f{LJ2)f{ri2,UJi,UJ2). (15) 

For hard-body interactions, this reduces to 

B2 = ^J duJ^duJ2f{^l)f{^2Wr\0l2). (16) 

where We"^(^i2) is the excluded volume between two rigid n-mers, which depends on the relative angle 612 between 
their axes. Analogous expressions hold for the monomer and dimer virial coefficients, li':^^ and B2^\ In particular, 

^(D) ^Pj^ f du;,du;2f{uJi)f{uj2)v^^\9i2), (17) 



(2) 

where p^) = nup is the number density of dimers and Ve (^12) is the corresponding excluded volume. Note that for 
rigid linear rods, the same angular distribution function, /(w), characterizes the full rod and any diatomic subunit of 
the rod. For the monomer fluid, 

5^) = ^.^, (18) 

where pm = n-Mp is the monomer number density. The monomer excluded volume vi^\ of course, has no angular 
dependence. Requiring that holds exactly then implies the following relation between the excluded volumes, 
independent of the angular distribution function and number densities, 

vi-Hdi2) = ai'MA'^ + ai^^nlvi^Hdu). (19) 

Using analytic results for the angle-dependent excluded volume between rigid n-mers first derived by Williamson 
and Jackson (WJ) |^ , we shall show that ( |l9| ) is indeed satisfied exactly for appropriate values of ai^^ and . The 
work of WJ was restricted to linear tangent hard-sphere chains, but it is straightforward to generalize their analysis 



to LFHS chains of arbitrary bondlength /. Details are contained in the appendix. The basic result is (see (A3)) 

vlr\e,2) = vi^Ho) + {n- iM\e^2), (20) 

where vi"\o) is the excluded volume for two parallel ri-mers {9i2 — 0), given by ( |Al| ) and (A2), and Vc'^\0i2) is 
the contribution from the so-called "central region" (in the terminology of WJ) of the overlap volume between two 
diatomic molecules. This is given by (A15) in the appendix. Note that, contrary to the remarks by WJ, the latter 
involves "compact and tractable" analytic formulae. 

Before comparing ( p^ ) and (pO|), one rewriting of the former is required. As already noted, the "reference" monomer 
and dimer fluids may be characterized by "effective" values of the numbers tim and no ■ Similarly, it may be allowed 
that the diameters and bondlengths of the reference particles differ from those of the original chain molecule. This 
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feature arises in GFD theory @ due to the assumption that the volume fractions of the reference fluids should equal 
that of the original n-mer fluid, rj = pvn, where w„ is the n-mer molecular volume. Equivalently, this requires that 
Vn = nMVi = nDV2, where vi and V2 are the molecular volumes of the reference monomers and dimers, respectively 
. The equality of the volume fractions of the reference fluids and n-mer fluid arises automatically in the case of 
tangent hard-sphere chains with I — d, and no adjustment of the number and sizes of reference particles is needed. 
This is no longer true when a chain of fused hard spheres (/ < d) is decomposed into reference monomers and dimers 
[0,^. Therefore, we suppose that the diameter of a reference monomer is denoted dM^ the diameter of a monomer in 
a reference dimer is denoted do, while the bondlength of the reference dimer is denoted In. It is these lengths which 
characterize the excluded volumes and We^''(0i2) in (|l9|), whereas all terms in ( pO[ ) for v'it^ (O12) are characterized 
by the diameter d and bondlength I of the original n-mer. Hence, it is appropriate to rewrite ( |l9| ) in terms of the 
original molecular dimensions. On dimensional grounds, as confirmed by the equations in the appendix, the dimer 
excluded volume scales as 



vP{012,dD,lD)^dl^ ( 012, 



do 



(21) 



where the dependence of v^e\0i2) on molecular lengths has been indicated. Then ( p9[ ) can be rewritten as 



(22) 



where We^-* and v\?\0i2) now refer to the original monomer and dimer inside the n-mer. (Note We^"* 
Comparing (pG) and (p2[) now gives 



47rd3/3.) 



{D) 



do 



(n - 1)^ 
no 



(23) 



and 



d Y b^^O) -(n- 1)2^(2) (0)] 



I'M 



(1) 
UMVe 



(2-3n)], 



(24) 



where the second line in (gj) follows from (Al) and (A2) in the appendix. Equations ( |23| ) and ( |24|) have been written 

in terms of the combined parameters Hmciti and riDak'^ because, as will be seen in the next subsection, only these 
combinations appear in the expressions for thermodynamic functions. As stated at the end of the previous subsection, 
and a'n^^ are seen to depend only on the geometric properties of the molecules. 
Of several prescriptions introduced in GFD theory [[7|j3l|] to adjust the effective parameter values, our reliance on 
the scaling relation (21) restricts us to those which conserve the reduced bondlength, lo/do = l/d = I*. These are 
the approaches analogous to the A, C, and AC versions of GFD theory [Q. However, it turns out that these methods 
all yield identical thermodynamic results in the present theory. This follows from the fact that, according to ( [2^ ) 

and (p4[), nMOii ^^ and nuai^"^ depend on the diameters and numbers (for fixed I*) only through the combinations 
nMid-M /d)^ and noidi) / dY , respectively, the values of which are uniquely determined by the condition of equal 
volume fractions 1^21. 



Let us briefly compare the present expressions for and al*^'' 



in (23) and ( p4[ ) with those given by GFD theory 
for LFHS chains. The latter were derived recently by Mehta and Honnell using GFD-A theory, which corresponds 
to taking dM = du = d while the efl^ective numbers of monomers Um and dimers no are adjusted to conserve volume 
fraction. For this case, Mehta and Honnell [|| showed that the function y„ in (^3|) has the value (n — 2) and thus (ITT 
reduces to 



(f) = -(n-2) 



= n - 1 {GFD - A) 



(25) 



independent of bondlength. For simplicity, we compare (|25| ) with ( |23| ) and (|24|) in the case of tangent spheres, for 
which Um — n and Ud — n/2. Eqs.(p3[) and (p3) then give 



5 



(M) (n-2)(ll-19n) 19 
^ W) ' 

n^al^) = ^±^-^2n, (26) 



n 



where we indicate the leading asymptotic dependence for large n. In this limit, the present parameters nj^io'n^^ and 
noaiF^ have nearly double the magnitude of those given by GFD theory. This leads, via to a significantly 



stronger angular dependence of the excluded volume (19) than would be predicted by GFD theory. One notes that 



al^^'' is negative for n > 2 according to both (^5|) and (|26D, and thus the non- ideal thermodynamic behavior of the 
fluid can be considered to involve a subtraction of the behavior of the monomer fluid from that of the dimer fluid. It 
can be shown pSf that the approximate factor of two difference between GFD and the present theory for the separate 

parameters aA/On and riDa'n''^ largely cancels out for the second virial coefficient in an isotropic phase, consistent 
with the findings of Mehta and Honnell |^ . By construction, the present theory yields exact values of B2 in both 
isotropic and ncmatic phases of LFHS chain fluids with arbitrary n. 



C. Free Energy Minimization 

The total Helmholtz free energy F for a uniform one-component LFHS chain fluid follows from (^, (^ij), and (pH): 
F ^ kTpV I dc./(c.)[ln(4^Ap/(c^)) - 1] + a^^'^ AF^'^'Xpm) + a^^) AF^^) (p,,; [/(c.)]). (27) 



The notation in (|2^) for AF^^^'^ and AF^^^ is meant to indicate that the former is a function of the monomer number 
density pM, while the latter is a function of pD while being a functional of the angular distribution f{uj). 

For the isotropic phase, we shall follow previous work on GFD theory and use the "exact" Carnahan-Starling |3^] 
and Tildesley-Streett ||5| equations of state to obtain Ai^(*^) and AF^^); 

^F(^Hpm) r?(4~3ry) _ 

Vkf " (1 - 7?)2 = PMacsiv), (28) 

= Puarsiv), (29) 

where 

F' = F + 3 = 4 + 0.37836r + 1.07860(r)3, 
G" = G - 3 = -2 + 1.30376r + 1.80010(r)^ 

H' = H -1^ 2.39803Z* + 0.35700(^*)^ (30) 

(and where F, G, H denote the same quantities as defined by Tildesley and Streett ||3^). As stated earlier, the volume 
fraction rj is assumed to be the same for the n-mer, dimer, and monomer fluids. In a uniform nematic fluid, the excess 
monomer free energy is still given by (p8|). The simplest conceivable ansatz for ls.F^^\pu', [/(w)]) is that suggested 
by decouphng theory p2| , p3|j29| ] , namely 

AF^^'^Pd: [f{^)]) = AF(^)(pz5)J[/(^)], (31) 

where J is given by 

From (|l7|), the functional J[f{u!)] is equivalent to the ratio between the second virial coefficients of the nematic and 
isotropic dimer fiuids. 

A self-consistent equation determining the equilibrium form of the angular distribution function follows by functional 
minimization of F with respect to f{uj), subject to the normalization condition J dLuf{u!) = 1. We obtain 
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Jduj'e{Lu')' 



where 



e{uji) — exp 



O (D) / N J duj2f{(^2We^\012) 

-2n£,a^ >aTs[m-r^ 7T\ 

JdL0,dL02vi''\di2)/W 



(33) 



(34) 



and we have used the fact that pD = nop- The excluded volume vi^\0i2) in ( |32D and (|3j) is that for reference dimers 
of diameter do and bondlength I n , but explicit factors of d'jj cancel in the ratio of integrals due to the scaling law 
( ^ ) and only the reduced bondlength I* is required in specifying those integrals. 

For a given volume fraction r], the integral equation (^3|) is solved numerically by iteration as in ref. |2^. To 
evaluate the angular integrals involving vi^\0i2) in (^J), we express the molecular orientation uji = (9i,4>i) in the 
director frame of reference, i.e., with polar angle 9i measured from the nematic director axis. Then /(wi) = f{Qi) is 
independent of the azimuthal angle (pi. The angle 9i2 between the two rod axes can be expressed as 



cos 6*12 = sin Oi sin 02 cos 02 + cos Oi cos 62 , 



(where we arbitrarily set =0). Then 



diJ2f{i^2)vi^H012)^ / sm02d92f{e2)v^^\Ou02), 
Jo 



(35) 



(36) 



where 



V^^H0i,e2) = / d^2vi^\9i2) ■ 



(37) 



Using the analytic formulae for vi^\9i2) from the appendix, the integrations over 62 and (j)2 in ( |36| ) and ( |37| ) are 
performed numerically by the trapezoid rule (for each pair of 61 and 62, the integration over 02 has to be done 
only once). The variables cos 6' and 4> are discretized on grids of typical stepsizcs A(cos6') = 0.005 and A(j> = 
7r/200. Changing the frame of reference to relative orientations (^12, 4'i2) could be applied to reduce the number of 
integration variables involved in the denominator of the exponential in ( P4| ) , but for numerical consistency and with 
a view towards cancellation of errors, that integral is evaluated here in a manner analogous to ( |36| ) (noting that the 
integration J c?wi/(47r) becomes redundant). Similar numerical integration techniques are used to evaluate J[/(a;)], 
the denominator of (33) and the order parameter S2 defined by (H) in Section [II A. 



Once the numerical solution of ( |33| ) is obtained, the corresponding equilibrium thermodynamic potentials can be 
evaluated from the Helmholtz free energy given by eqs.(p7|)-(3^). Substituting ( |3^ ) into the logarithm of (27), one 
obtains 



F 
kTV 



= P 



In 



AnpA 
J dhjeiuj 



- 1 



The chemical potential p, follows from 



Hence 



kT 



In 



A-Kpk 



J dwe{uj) 



nMa\ ' — — — +nDa\, V[/(w)] 



A/ig(??) 
kT 



2aTs{v) 



(38) 



(39) 



(40) 



where ApM and Apu are the excess chemical potentials of the monomer and (isotropic) dimer fluids, respectively: 

A/iM(?7) 



kT 
Amd(??) 
kT 



acsiv) +Va.'csiv), 
arsiv) +m'Tsiv), 



(41) 
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where a'(?7) = da{r])/dri. Finally, the equilibrium pressure is obtained from P = pii — F/V, which gives 

^=P + PMai"''>AZM{v) + PDai''^J[f{u;)]AZD{v), (42) 

where AZm and AZo are the excess monomer and (isotropic) dimer compressibility factors, respectively, given by 

2ry(2-77) 



azm{v) = m'csii) 



3 



, , . rj{F' + G'ii- H't?) , , 

AZu{ri) = ^o:^.s(r|) = ^ ^^_;^^3 (43) 

The isotropic limits of the preceding equations are obtained by setting /(w) = l/(47r). In this limit, J\J(uj)\ = 1. 
Coexistence between isotropic(I) and nematic(N) phases is then evaluated by solving the equations P{rii) — P{tin) 
and fJ.{rii) = fiiji^), which is done by a Newton-Raphson procedure. 

III. RESULTS 
A. LTHS Chains 

We consider first the case of tangent hard-sphere n-mers, before examining n-mers with arbitrary bondlength. The 
primary limitation to our analysis of the theory is the lack of extensive simulation data available for comparison. Of 
the Monte Carlo studies which have examined hard-sphere chains, the majority have focused on semi-flexible systems 
of tangent chains. Several of these studies have dealt with the limit of infinite rigidity and thereforeyield results 
which can be used here for comparison. In particular, the recent studies by Yethiraj and Fynewever ||ll|,^ provide 
us with substantial analyses of the 8-mer and 20-mer LTHS chains. Williamson and Jackson recently examined 
7-mer LTHS chains and performed comparisons with theoretical treatments involving the exact excluded volume for 
the chains | |25[ |. These studies will be the primary sources for comparisons with our results. 

The conventional orientational order parameter 5*2 is defined as 

S2 = J P2{cos 9) f{Lo)diu, (44) 

where P2 is the second Legendre polynomial. This parameter has the limits of 52 = for an isotropic phase and 
^2 = 1 for a perfectly aligned nematic phase. Figure 1 (a) and (b) show the order parameter ^2 as a function of 
volume fraction 77 for 8-mers and 20-mers, respectively. The Monte Carlo simulation data is obtained using both 
constant-pressure (NPT) and constant-volume (NVT) methods, which are seen to be in very good agreement. This 
data is compared with the present theory and with the Parsons theory p^ ], as employed by Yethiraj and Fynewever 
pT| , pl| (using a simple hard-sphere reference fluid and incorporating the excluded volume derived by Williamson and 
Jackson p5|]). 

In general, the present theory is in very good agreement with the simulation data for the degree of ordering in each 
system. The order parameter at which the nematic phase flrst appears is of the order ^2 > 0.5. The degree of ordering 
continues to increase dramatically over an extremely short density range, until ^2 « 0.8, after which the value of the 
order parameter begins to level off. This general behavior is also evident in the Parsons theory, although the values 
obtained for S2 and the coexistence densities are less accurate. For the 8-mer LTHS chain system, the values of the 
order parameter predicted by our theory are in excellent agreement with the values obtained by simulation. Our 
theoretical results for the 20-mers exhibit slightly lower values of 5*2 in comparison with the simulations, although the 
amount of data available is much smaller in this case. The Parsons theory systematically underestimates the degree 
of ordering in the nematic phase, as is evidenced in Figs. 1(a) and (b). 

In Figs. 2(a) and (b), the reduced pressure P* — Pvn/kT is plotted as a function of volume fraction for the 8- 
mer and 20-mer LTHS chains, respectively. The simulation data clearly show the coexistence densities between the 
isotropic and nematic branches. The values for these densities as yielded by both theory and simulation are given in 
Table 1, along with the order parameter S2 and P* in the coexisting nematic phase. It is seen that the coexistence 
regions predicted by the present theory agree very well with simulation, although the reduced pressure at any density 
shows poor agreement. The present theory overestimates the value of P* throughout the density range, a discrepancy 
which becomes more pronounced at higher densities. On the other hand, the Parsons theory predicts coexistence 
densities substantially higher than those observed in simulation, as well as a larger coexistence range, as is depicted 
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on both graphs. The Parsons theory underestimates the values of the reduced pressure as a function of density for 
both systems. Yethiraj and Fynewever note that SPT treatments also underestimate the pressure for the 8-mer 
system, while overestimating it for larger molecules. 

The 7-mer LTHS chain has been recently examined by Williamson and Jackson [WJ] ||2^] in light of the explicit 
calculation of the excluded volume for LTHS chains WJ performed Monte Carlo simulations for a system of 

N = 576 molecules. The results of the simulation were compared with several theories which were modified to 
incorporate the calculation of the exact excluded volume. In particular, WJ compare the simulation results with 
those obtained from a modified version of the Vega and Lago theory The original Vega and Lago theory is 

based on the form of the decoupling approximation in (|^) and (^), where the reference fluid is the isotropic phase 
of the n-mer fluid. WJ modify the original Vega and Lago theory to incorporate their exact analytic calculation of 
the excluded volume into B2, and obtain the isotropic phase free energy from TFT [|l]j36). In Fig. 3, we compare 
the results for S2 from the present theory with the Monte Carlo simulation data. As in Fig. 1, the present theory 
predicts the order parameter with excellent accuracy. The general behavior of S2 is similar to the previous graphs. 
The comparison of reduced pressure versus volume fraction for the 7-mer system is made in Fig. 4. The modifled 
Vega and Lago theory yields accurate pressures in comparison with simulation, but the predicted transition densities 
(evidenced by the plateau in the trace) are too low. This is also indicated in Table 1, which shows that the present 
theory yields coexistence densities lying within the ranges obtained through the simulations. The Vega-Lago theory 
also overestimates the value of iS'2 in the nematic phase. However, as in Fig. 2, the values of reduced pressure given 
by the present theory clearly exceed the rep orted simulation data. This discrepancy appears to be systematic in the 
theory and will be discussed in Section It should be noted here that the WJ simulation data predict a smectic 
phase to occur at volume fractions greater than 0.37. 

The relationship between the length of molecules and the range of the coexistence region is shown in Fig. 5. The 
general trend is an increasing difference between the densities of the isotropic and nematic branches at coexistence as 
the number of monomers constituting the molecule increases. This trend is consistent with simulation, as is shown in 
this figure for 7-mers, 8-mers and 20-mers. The Parsons theory and various SPT treatments have yielded a similar 
trend [0, although these theories tend to overestimate the width of the coexistence region. In addition, the values 
of the coexistence densities obtained from these theories are in poor agreement with simulation, while the present 
theory provides a quantitatively accurate description of these densities. 



B. LFHS Chains 

Simulations of LFHS chains have been performed by Whittle and Masters for the cases of 6-mers with a reduced 
bondlcngth of I* = 0.5 and 8-mers with I* — 0.5 and 0.6. These three systems were recently investigated by Mehta 
and Honnell 1^ using GFD theory, comparing their results with the Whittle and Masters simulation data and 
with a modification of TPT 0]. The GFD theory and the modified TPT are unable to treat ordering, and therefore 
are inappropriate at densities above the isotropic-nematic transition which was evidenced in the simulation of the 
8-mer system with I* — 0.6. In Fig. 6, we perform an analogous study of the three fluids, comparing the previous 
results [|j with those of the present theory. Figure 6(a) shows the reduced pressure as a function of volume fraction 
for the 6-mer LFHS chains with I* = 0.5. There is no evidence of a nematic transition in the simulations, nor is there 
any indication through the theoretical treatments. The present theory and the modified TPT both overestimate the 
values of the reduced pressure, while the GFD theory appears to predict the reduced pressure to a fair degree of 
accuracy. This trend is also apparent in Figs. 6(b) and (c) for the 8-mer cases. In the latter plots, the present theory 
is seen to predict a stable nematic branch at sufficiently high volume fractions. The simulation data in Fig. 6(b) 
exhibit no such transition, and we are unaware of any other simulations which have been done for this system. It 
should be noted that the monomeric reference fluid does not exist as a stable fluid above rj « 0.494 ||] and as such 
the theory is probably invalid at these high densities. 

Figure 6(c) considers the 8-mcr fluid with I* — 0.6, a system which clearly exhibits nematic ordering. At first glance 
it appears that the present theory, while succeeding in predicting a stable nematic branch, predicts that it occurs 
at volume fractions much greater than the simulations indicate. However, it should be noted that although Whittle 
and Masters \ ^7\ report a nematic branching at f« 0.33, the order parameter was not found to be stable until much 
higher volume fractions, rj > 0.419. The predictions of the present theory for transition properties are compared with 
the simulation results in Table 1. Once again, we see that there is good agreement for the value of the nematic order 
parameter 6*2. 

It is of interest to examine the present theory in the "spherocylinder" limit, which is obtained in the limits n — > 00, 
I —^ 0, such that L = [n — 1)1 remains finite. As discussed in the appendix, our analytical result for the LFHS excluded 
volume reduces (as should be expected) to that for hard spherocylinders of cylinder length L in this limit. Figure 7 
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shows the nematic and isotropic coexistence volume fractions for LFHS n-mers as a function of n, for a fixed value of 
L/d — 19, in comparison with the values for spherocylinders obtained from the Lee theory and from Monte Carlo 
calculations (using the Gibbs-Duhem integration procedure) by Bolhuis and Frenkel [Q. One sees that both the 
individual volume fractions and their differences increase slowly with n, and approach the asymptotic spherocylinder 
values reasonably well, although underestimating the coexistence width in the limit. In more detail, it can be shown 
that the present theory in the spherocylinder limit yields exactly the same angular distribution function f{uj), at a 
given volume fraction 77, as in the Lee theory ||2^ . However, the free energy and pressure are actually predicted to 
diverge in this limit, which accounts for the narrower coexistence width found in Fig. 7. 



IV. CONCLUDING REMARKS 



In this paper, we have introduced two key modifications into the generalized Flory-dimer (GFD) theory to describe 
nematic behavior in hard-sphere chain fluids of arbitrary intramolecular bondlength. The first modification is the 
inclusion of the exact excluded volume and second virial coefficient for LFHS chain molecules, based on generalizing 
the earlier calculation of Williamson and Jackson [^5| . This procedure results in mixing parameters for the reference 
monomer and dimer fluids which are dependent only on the geometric properties of the molecules (see (^3|) and 
(p^)). A related feature is that the ansatz of equal volume fractions uniquely determines the values of the relevant 
combinations of the effective reference parameters nu,du,nM and dM, in contrast to previous studies based on GFD 
and related theories Pj7|,|3l|]. The second modification of GFD theory, in order to account for the possibility of 
nematic ordering in the system, is the weighting of the excess dimer free energy by the ratio between the second virial 
coefiicients of the nematic and isotropic dimer fluids (see ( ^ and (|3|)). This ansatz, suggested by decoupling theory, 
is dependent on the angular distribution function, which is determined self-consistently. 

The present theory is found to be in excellent agreement with simulations of LTHS fluids in determining coexistence 
densities and the nematic order parameter as a function of density. As is seen clearly in Table 1 and Fig. 5, the 
transition densities for the 7-mer, 8-mer and 20-mer fluids predicted by the theory fall within the coexistence range 
found in simulations. In particular, the 7-mer system studied by Williamson and Jackson ||20|| shows much better 
agreement with the present theory than with the Vega and Lago |l^ and the Parsons theories for the coexistence 
densities as well as the nematic order parameter at the transition. The present theory is able to account for stable 
nematic branches in LFHS chains of sufficient length, for which theoretical work has been lacking up to now. The 
8-mer with I* = 0.6 simulated by Whittle and Masters yielded transitions at lower densities than predicted by the 
present theory, although the simulated nematic branch did not become stable until densities similar to those predicted 
in our work. In addition, transitions between the isotropic and nematic phases are predicted by the present theory 
for several systems which have yet to be simulated (see Figs. 5 and 7). The agreement between the present theory 
and simulations is promising and it is hoped that further work will be encouraged. 

For all systems studied, the present theory yields values for the pressure which exceed those given by available 
simulation data as well as previous theories. This discrepancy is no t acc ounted for within this paper, and further 
work needs to be done to resolve the problem. It was noted in Section III B that the free energy and pressure given by 
the present theory actually diverge in the spherocylinder limit n —^ 00, I — *■ 0, with (rt — 1)^ = L finite. This divergence 
can be traced to the incorrect limiting behavior of the Tildesley-Streett (TS) dimer equation of state for small 
reduced bondlength I*. The TS equation of state predicts that the reduced isotropic second virial coefficient (as 
well as higher virial coefficients) varies linearly with I* (see (|30|)), in contrast with the correct leading-order variation 
proportional to (l*)^ (see (A17)). However, use of other dimer equations of state which do have the correct dependence 
for I* — > 0, such as the "improved scale particle theory" of Boublik and Nezbeda produces negligible differences 
from most of the results described in this paper, apart from those shown in Fig. 7 for n > 40. 

Despite the mixed agreement of our results with those of available simulations, we believe that the approach 
described here has advantages over other current density functional theories of chain fiuids [p,p yTl| , p^ -pT[|, particularly 
for c onsidering extensions of the theory to non-uniform fluids and ones containing semi-flexible molecules. In Section 



11 A , we briefly indicated these generalizations of the theory, which we are currently investigating in more detail. 
Our approach retains the geometrically motivated spirit of the GFD and TPT theories in utilizing reference fluids 
composed of monomer and dimer subunits, whose properties can be readily determined. This contrasts with the form 
of the decoupling approximation used in refs. [ p"9|j2C| ], requiring a priori information about the thermodynamics of the 
isotropic phase of the full system being considered, which may be either unavailable or computationally difficult to 
obtain in the more general cases. At the same time, the dimer reference fluid has the crucial property of involving 
orientational degrees of freedom. This feature corrects a limitation of previous density-functional approaches for 
non-uniform chain fiuids P,p|Jl^-p^ based on a purely monomeric reference fluid, which are unable to account for 
orientational ordering effects. 
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APPENDIX : EXCLUDED VOLUME OF LFHS CHAINS 



Williamson and Jackson (WJ) |25j recently derived an expression for the excluded volume between a pair of LTHS 
n-mers at an arbitrary relative orientation 612. Here we extend their analysis to LFHS n-mers of arbitrary bondlength. 

The excluded volume for parallel orientation, 9i2 = 0, is a straightforward extension of the WJ result, and is 
described by a chain of (2n — 1) overlapping spheres of radius d and volume vs — 4,71(1^/3 separated by the rt-mer 
bondlength I. The volume is given by the relation 

(0) = {2n - l)vs - 2{n - l)vo. (Al) 

The overlap volume between a pair of adjacent spheres, Vo, is dependent on the bondlength-to-diameter ratio l/d, and 
is given by 



,3 

.0 = 



3 / A 1 / r ^ 



(A2) 



For a pair of 71-mer molecules in an arbitrary orientation 612, the excluded volume is represented diagrammatically 
as overlapping spheres of radius d whose centers lie on the corners of a rhombus (Fig. 8). The centers of adjacent 
spheres are separated by I and the rhombus has angle 6*12. We define the "central region" (in the terminology of WJ) 
as the parallelpiped based on the rhombus, taken to lie in the xy plane, and extending along the z axis to distances 
±d. As in WJ |2^, it can be shown that the excluded volume outside the central region is equal to the excluded 
volume in the 9i2 — orientation. Furthermore, one can show that the excluded volume contained within the central 
region for an n-mer is {n — 1)^ times that of a dimer. The total excluded volume is thus given by 



vi'^^ (012) = vi''^ (0) + {n- ifvP {0,2), (A3) 



where Ve"' (0) is given by ( Al) and Vc^' (^12) is the contribution from the central region to the excluded volume between 



„(2) 



two diatomic molecules. 

Following WJ |2^, it is convenient to evaluate vi^\9i2) by considering infinitely thin slices parallel to xy at varying 
heights z. Due to mirror symmetry in the xy plane, we need to examine only distances < z < d. There are three 
distinct ranges of z which must be considered. The outermost range of z is characterized by the absence of two-body 
overlaps between the circular cross-sections through the spheres. The total excluded area of a slice parallel to xy 
within this range is that of a circle, 

Ai{z)^n{d'~z'), (A4) 

where z ranges from d to /i//, the height at which two-body overlaps begin to occur. The distance hn is dependent on 
the magnitude of the angle between the molecules, 612- For large enough 6*12 (7r/3 < 6*12 < 7r/2), overlap between the 
circular cross-sections first occurs for adjacent pairs of spheres. At this height, the radius of a circular cross-section 
is r = 1/2, and the corresponding value of 2; = hn is: 

4y = ^Jd^-{l/2f. (A5) 

For smaller 612 (0 < 9i2 < 7r/3), the first two-body overlaps occur between opposing spheres (Fig. 9). This height is 
given by the relation 



hf^ = ^d^-Psm^0i2/2). (A6) 

For both angular cases evaluated above, the region of two-body overlaps extends until z = hm, at which height the 
circular cross-sections through the spheres begin overlapping in a three-body configuration. The value for this height 
is given by 
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T-III 



I 



^2 COS (012/2)^ 

Throughout the range hm < z < d, the total excluded area of a parallel slice can be expressed as 



(A7) 



(A8) 



where A^^{z) is given by (A4) while Ao{z) is the total two-body overlap area between the circles at height z. This 
overlap area has distinct contributions depending on whether the overlapping spheres are adjacent or opposite from 
each other (Fig. 9), and is given by 



Ao{z) = 2ao{l, z) + ao{l, z), 



(A9) 



where / = 2Z sin (6*12/2). Here, ao{l' , z) is the area of overlap between a pair of circles of radius r, separated by 
given by 



aoil',z) = 2 



9 —1 
r cos 



;/2\ P 



(AlO) 



keeping in mind that r — y/cP 



The relation ( AlO ) differs from equation (7) of WJ by a factor of 2 because 



WJ calculate only half of the overlap area. Remembering that the contributions of the two types of overlap begin at 
different values of z, the total contribution to the central excluded volume from regions / and // is 



dz A[' (z) 



id-hiii)--{d^~h 



hi) 



dz ao{l, z) 



dz ao{l, z) 



(All) 



where the additional factor of 2 in (All) is due to mirror symmetry. 

tract 

l'/2 



The integrals involving a'^s are quoted by WJ as being "intractable" , but we have derived a compact and tractable 
form for the indefinite integral: 



dz' ao{l',z') = 2{d^z- z^/3) 



cos 



d^-z^- {i'/2y 



1/2 



4' 



V^'-(^'/2)' 



'-d^ tan^i 



Vz/{2d) 



rf2 



{V/2Y 



1/2 



(A12) 



One can verify that the z-derivative of the right-hand side of ( A12| ) equals ao{l' , z). 

The innermost part of the central region is defined by z < hm, where three-body overlaps occur. The whole volume 
of this region is excluded, and thus its contribution from this region to Vc is equal to 



vi''{0i2)^2l'd 



sm 612 



sin^ (012/2) 



1/2 



The total excluded volume of the central region is then 

vf\e^2) = v'/\e,2) + vi'\e^2), 



(A13) 



(AM) 



where the various contributions are given by (All) to (A13) for ^12 < 7r/2. Due to chain inversion symmetry, the 



excluded volume for ^12 > 7r/2 is given by Vc{0i2) — Vc'^\'k — 612). Substituting the expressions for hYj,hYj and 



(1) 1,(2) 
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hill from (A5), (|A6|) and (|A7j ) into (|All| ) and ( [A12| ), and applying several trigonometric identities (most importantly 

the relation tan~^(yl) — tan^^(i?) — ta.'Dr^[{A — B)/{1 + AB)]), we are able to reduce the result for v'"c\9i2) to the 
following compact form: 



v'i'\0^2)^l^f^d^ sin (012/2) 5(012) 



+4^ d^^'-(l 
Y 3 V2 

+4/ sin (6*12/2 
-d^ 1 tan"^ 



tan 



I sin (6I12/2) 



d2 _ ^ sin^ (012/2) 



tan 



I cos(0i2) 

5(^12) 



sin (012/2) 5(012) 
4(i2-P(l + 2sin2 (012/2)) 



where 



5(01 



4cos2 (012/2) 



d2 



1/2 



(Al5a) 



(A15b) 



It can be verified that the preceding result for the angle-dependent excluded volume of LFHS n-mers reduces to 
that for hard spherocylinders in the li mit t hat n — » c x3, I ~* 0, and L = [n ~ \)l remains finite, where L is identified 
with the cylinder length. In this limit, (Al) and (A2) give 



,(") 



(0) 



47rd3 



^+2d 



where vsc is the spherocylinder molecular volume. From ( Al5a| ,b), one finds readily that 

The total excluded volume ( |A3| ) then becomes 

v'i'\ei2) 8wsc + 2L2dsin0i2, 
which agrees with the result derived for hard spherocylinders p3|,p4|. 



(A16) 



(A17) 



(A18) 
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TABLE 1. Coexistence results from simulation and theory. 
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FIG. 1. Variation of order parameter 52 with volume fraction r/, comparing present theory with the Parsons theory and 
Monte Carlo data[ll,21], for (a)8-mer LTHS chains, and (b)20-mer LTHS chains. 



FIG. 2. Variation of the reduced pressure with volume fraction 77, comparing present theory with the Parsons theory and 
Monte Carlo data[ll,21], for (a)8-mer LTHS chains, and (b)20-mer LTHS chains. 



FIG. 3. Order parameter 52 vs. volume fraction for LTHS 7-mers, comparing present theory with Monte Carlo data[20]. 

FIG. 4. Reduced pressure for LTHS 7-mers as a function of volume fraction rj, comparing present theory with Monte Carlo 
data and the modified Vega and Lago theory from Ref. [20] . 

FIG. 5. Comparison of volume fractions of the isotropic and nematic phases at coexistence for LTHS n-mers, as a function 
of the number of monomers n, between the present theory and Monte Carlo simulation data. 

FIG. 6. Comparison of the reduced pressure between this work, the Mehta and Honnell GFD theory and the modified TPT 
in Ref. [8], and Monte Carlo simulations[37] for (a)LFHS 6-mers with bondlength to diameter ratio 1* = 0.5, (b)LFHS 8-mers, 
r = 0.5, and (c)LFHS 8-mers, l* = 0.6. 

FIG. 7. Volume fractions of the isotropic and nematic phases at coexistence as a function of number of monomers n for 
LFHS chains of constant length with L/d — 19. The spherocylinder limits are given by the Lee theory [23] and by the Monte 
Carlo simulation results of Bolhuis and Frenkel[38]. 



FIG. 8. Diagrammatic representation of the excluded volume for n = 2. The slice is taken through the z — plane, where 
the radius of each circle is the monomer diameter, d. 



FIG. 9. Second case for two-body overlap, when < 612 < 7r/3. 
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Figure 1(a). K.M. Jaffer: J. Chem. Phys. 
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Figure 1(b). K.M. Jaffer: J. Chem. Phys. 
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Figure 2(a). K.M. Jaffer: J. Chem. Phys. 
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Figure 2(b). K.M. Jaffer: J. Chem. Phys. 
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Figure 3. K.M. Jaffer: J. Chem. Phys. 
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Figure 4. K.M. Jaffer: J. Chem. Phys. 
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Figure 5. K.M. Jaffer: J. Chem. Phys. 
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Figure 6(a). K.M. Jaffer: J. Chem. Phys. 




Figure 6(b). K.M. Jaffer: J. Chem. Phys. 
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Figure 6(c). K.M. Jaffer: J. Chem. Phys. 
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Figure 7. K.M. Jaffer: J. Chem. Phys. 
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Figure 8. K.M. Jaffer: J. Chem. Phys. 




Figure 9. K.M. Jaffer: J. Chem. Phys 




